A FAMILY OF MARTINGALES GENERATED BY 
A PROCESS WITH INDEPENDENT INCREMENTS 



Josep Llui's Sole and Frederic Utzet 

Abstract. An explicit procedure to construct a family of martingales generated by a process with independent 
increments is presented. The main tools are the polynomials that give the relationship between the moments 
and cumulants, and a set of martingales related to the jumps of the process called Teugels martingales. 



1. Introduction 

In this work we present an explicit procedure to generate a family of martingales from a process 
X = {X t , t > 0} with independent increments and continuous in probability. We extend our results 
exposed in [8], where we dealt with Levy processes (independent and stationary increments); in that 
case, the martingales obtained were of the form M t = P{X t , t), where P(x, i) is a polynomial in x and i, 
and then they are time-space harmonic polynomials relative to X . Here, the martingales constructed are 
polynomials on Xt but, in general, not in t. Part of the paper is devoted to define the Teugels martingales 
of a process with independent increments; such martingales, introduced by Nualart and Schoutens [5] for 
Levy processes, are a building block of the stochastic calculus with that type of processes. 

2. Independent increment processes and their Teugels martingales 

Let X = {X t , t > 0} be a process with independent increments, X Q = 0, continuous in probability 
and cadlag; such processes are also called additive processes, and we will indistinctly use both names. 
Moreover, assume that X t is centered, and has moment of all orders. It is well known that the law of 
Xt is infinitely divisible for all t > 0. Let of the variance of the gaussian part of Xt, and Vt its Levy 
measure; for all these notions we refer to Sato [6] or Skorohod [7]. 

Denote by v the (unique) measure on 6((0, oo) x Ro) defined by 

u((0,t] x B) = v t {B), BeB{R Q ), 

where Ro = R — {0}. By an standard approximation argument, we have that for a measurable function 
/ : R — > R, and for every t > 0, 

|/(x)| v(ds, dx) < oo / \f{x)\ v t {dx) < oo, 

(0,t]xK o ilo 



and in this case, 



f{x)v{ds,dx)— j f(x)vt{dx). 

(0,t]xR 
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Note that since for every t > 0, v t is a Levy measure, then v is cr-finite. To prove this, observe that 
^t({|x| > 1}) < oo, and v t ((l/(n + 1), 1/n]) < oo and v t ([—l/n,—l/(n+l)))<oo,n>l. So there is 
a numerable partition of Ro with sets of finite v t measure, Vt > 0. Then, we can construct a numerable 
partition of (0, oo) x R , eacn se t with finite i7-measure. 

Write 

N(C) = #{t : (t, AX t ) G C}, C G B((0, oo) x Ro), 

the jump measure of the process, where AX t = X t — X t - ■ It is a Poisson random measure on (0, oo) x R 
with intensity measure v (Sato [6, Theorem 19.2]). Define the compensated jump measure 

dN(t, x) = dN(t, x) - du(t, x). 
The process admits the Levy— Ito representation 



X t = G t + // xdN(t,x), (2) 

J J (0,t] xR 

where {G t , t > 0} is a centered continuous Gaussian process with independent increments and variance 

It is also well known the relationship between the moments of an infinitely divisible law and the 
moments of its Levy measure (see Sato [6, Theorem 25.4]). In our case, as the process has moments of all 
orders, for alH > 0, 

/ \x\ v t {dx) < oo and / \x\ n v t {dx) < oo, Vn > 2. 

J{|a;|>l} JR 

Write 

F 2 (t) = cr t 2 +/ x 2 v t (dx) and F n (t) = x n v t {dx), n > 3. (1) 
Since J^y^ \x\ vt{dx) < oo and E[X t ] = 0, the characteristic function of X t can be written as 

<t> t (u) = exp { - - a\ u 2 + I (e lux - 1 - iux) v t (dx) } . 
2 Jr 



It is deduced that for n > 2, F n (t) is the cumulant of order n of X t (for n = 2, E[Xf] = F 2 (t)). Also of 
is continuous and increasing (Sato [6, Theorem 9.8]). 

Proposition 1. The functions F n (t), n>2, are continuous and have finite variation on finite intervals, 
and for n even, they are increasing. 

Proof. 

Consider < u < t < v, and write U = [u,v]. From the continuity in probability of X, 

lim X™ = X™, in probability. 

s—>t,sEU 

Moreover, Vs G U, \X S \ < sup re(7 \X r \, and since A" is a martingale, by Doob's inequality, 

E[sup \X r \ n ] < CsupE[|X r |™] < CE[\X v \ n ] < oo. 

reu reu 
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So by dominated convergence it follows that the function 1 1— > E[X™] is continuous. Since the cumulants 
are polynomials of the moments, it is deduced the continuity of all functions F n (t). 

To prove that F n (t) has finite variation on finite intervals, consider a partition of [0, t]: < to < ■ ■ ■ < 
t k = t. Then 

V \F n {tj) - F n (tj-i)\ =J2\ x n v(ds,dx)\ 

j=l j=i J -/(tj-i,tj]xR 

<V// \x\ n v(da,dx) = ff \x\ n v(ds,dx) < oo. ■ 



Consider the variations of the process X (see Meyer [4]): 

xW = [X,X] t = *>+ £ (M s ) 2 

0<s<t 

= £ (A*.)", n>3. 

0<s<t 

By Kyprianou [2, Theorem 2.7], for n > 3 (the case n = 2 is similar), the characteristic function of X^ 
is 

exp { ff (e luxn - 1) u(ds, dx) } = exp { f (e mx - l) v\ n) (dx) } , 

where v[ n ^ is the measure image of v t by the function x ^> x n , which is a Levy measure. So X^ has 
independent increments. Also by Kyprianou [2, Theorem 2.7], for n > 2, 

E[Af l) ] = F n (t) and E[(X t (n) ) 2 ] = F 2n (i) + (F„(t)) 2 . 

Therefore, combining the independence of the increments and the continuity of F n (t) , it is deduced that 
X^ is continuous in probability. 

By Proposition 1, F n (t) has finite variation on finite intervals. Hence, the process 

X^=F n (t) + (X^-F n (t)) 

is a semimartingale. 



The Teugels martingales introduced by Nualart and Schoutens [5] for Levy processes can be extended 
to additive processes. In the same way as in [5], these martingales are obtained centering the processes 

Y t V=X t , 

Y } n) = JTW - F n (t), n>2, 
They are square integrable martingales with optional quadratic covariation 

[y(") j y( m )] t = x'™ +m \ 
and, since F 2 „(i) is increasing, the predictable quadratic variation of Y^ is 

(Y^) t =F 2n (t). 
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3. The polynomials of cumulants 

The formal expression 

exp^ ^>„— | = 2>„ — . (3) 

^ n=l ' n=0 

relates the sequences of numbers {n ni n > 1} and {/x„, n > 0}. When we consider a random variable 
Z with moment generating function in some open interval containing 0, then both series converge in a 
neighborhood of 0, and (3) is the relationship between the moment generating function, ip(u) = E[e" z ], 
and the cumulant generating function, logip(u). Moreover, /i„ (respectively, n n ) is the moment (respec- 
tively, the cumulant) of order n of Z, and the well known relationship between moments and cumulants 
can be deduced from (3). The first three are 

Mi = «i, 

M2 = re? + re 2 , 

= ref + 3«ik 2 + k 3 , . . . 

If the random variable Z has only finite moments up to order n, the corresponding relationship is true 
up to this order. 

There is a general explicit expression of the moments in terms of cumulants in Kendall and Stuart 
[1], or formulas involving the partitions of a set, see McCullagh [3]. In general, \i n is a polynomial of 
Ki,. .. ,k„, called Kendall polynomial. Denote by r„(xi, . . . , x n ), n > 1, this polynomial, that is, we have 

[i n = r„(Ki, . . . , K n ). 

Also write r = 1. These polynomials enjoy very interesting properties, as the recurrence formula that 
follows from Stanley [9, Proposition 5.1.7]: 

r„+i(ari, . . .,x n+ i) = Y ( • ) T j( x i' ■ ■ -,Xj)x n+1 -j. (4) 
3=0 ^ 

We also have 

^ = y^.jT n _j(xi, . . . ,x n -j), j = l,...,n. (5) 

Computing the Taylor expansion of T n (x\ + y, X2, ■ ■ ■ , x n ) at y — 0, we get the following expression that 
we will need later: 



T„(xi +y,x 2 ,. ,x n ) = ^ ( . JT n - j (x 1 , . . ..x^^y . (6) 
Interchanging the role of x\ and y, and evaluating the function at we obtain 



r„(xi,a;2, ...,x n )=^2 f • ) r n-j(0, 

3=0 



X2,...,X n -j)x{. (7) 
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4. A FAMILY OF MARTINGALES RELATIVE TO THE ADDITIVE PROCESS 

The main result of the paper is the following Theorem: 
Theorem 1. Let X be a centered additive process with finite moments of all orders. Then the process 



M t (n) = r, 



,(x t ,-F 2 (t),...,-F n (tj) 



is a martingale. 
Proof. 

Let n > 2. We apply the multidimensional Ito formula to the semimartingales X t , F 2 (t), . . . ,F n (t). 
By Proposition 1, the functions F 2 (t), . . . , F n (t) and of are continuous and of finite variation. From (5) 
and the fact that [X, X]^ — of and [Fj, Fj]^ = 0, we have 

M<"> = nj* Mir 1} dX s g (") f Q M^dF 3 {s) + \n(n 1) J* d^) 

+ (T n (X s _+AX s ,-F 2 (s),...,-F n (s))-T n (X s _,-F 2 (s),...,-F n ( S )) 

0<s<t 

- nAX s r n _! (X s _, -F 2 (s), . . . , -F n (s)) 

Applying (6), 



T n (X s _ + AX S , -F 2 (s), -F n (s)) = J2 QM^- j) (AX s ) j . 

(n) 

Then, the jumps part given in the expression of M t v is 

e e ( ■) m ^ 3) ( ax *y = e ( ■) jf M s-~ j) - (;) jf* Mi—) ^ 



J=2 

Therefore, 



Mi"-2) d(<T 2) 



M^=p^J*M^_- j) dY^. (8) 

Moreover, (Mi*-*) 2 is a polynomial in X t , F 2 (t), . . . , Fk(t). Taking expectations and using the relationship 
between moments and cumulants, and that the cumulants of X t are F n (t), n > 2, we obtain that 

E[(M^) 2 } =P(F 2 (t),...,F 2k (t)), 
for a suitable polynomial P. Then, for every t > 0, we have 

E[ f (M { s k J) 2 d(Y^) s ] = f ' E\(M^) 2 ]dF 2j (s) = f P(F 2 (s),...,F 2k (s))dF 2j (s) <oo, 
Jo Jo 1 J Jo 

So all the stochastic integrals on the right hand side of (8) are martingales. ■ 
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Remark 1. It is worth to note that it follows from the preceding Theorem that the function 

g n (x,t)=T n (x,-F 2 (t),...,-F n (t)) 
is a time-space harmonic function with respect to X t . By ( 7) 

n 

g n (x, t) =^n-j(0, -F 2 (t), F n (t))x j . 

3=0 

In general, g n {x,t) is a polynomial in x. If F n (t), n > 2, are polynomials in t, then g n {x,t) is a time- 
space harmonic polynomial; this happens for all Levy processes with moments of all orders, and for some 
additive process; see the example below. 

Example. Let A(t) : M+ — ► M+ be a continuous increasing function, and J be a Poisson random 
measure on R + with intensity measure fi{A) = J A A(dt), A G B(R+). Then the process X — {X t , t > 0} 
defined pathwise, 

X t (w)= [ J(ds,w)-A(t) 
Jo 

is an additive process; it is a Cox process with deterministic hazard function A(t). From the characteristic 
function of X t we deduce that the Levy measure is 

v t (dx) = A(t)Si(dx), 

where Si is a Dirac delta measure concentrated in the point 1. Hence, 

F n (t)=A(t), n>2. 

Note that the conditions that we have assumed on A are necessary to obtain an additive process, but it 
is not necessary (though not very restrictive) to assume that A is absolutely continuous with respect to 
the Lebesgue measure. 

The function defined in Remark 1 is 

g n (x,t) =T n (x,-A(t),...,-A(t)). 

Hence, when A(t) is a polynomial, g n {x,t) is a time-space harmonic polynomial. 

Denote by C n (x,t) the Charlier polynomial with leading coefficient equal to 1, then (see [8]) 

n 

g n (x,t) = J2^ n) Cj(x,A(t)), 

3=1 

where = 1 and 

Aa = E(")A« * = l,...,n-l. 

3=k V/ 
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